Geometric phase for open quantum systems and stochastic unravelHngs 
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We analyze the geometric phase for an open quantum system when computed by resorting to 
a stochastic unravelUng of the reduced density matrix (quantum jump approach or stochastic 
Schrodinger equations). We show that the resulting phase strongly depends on the type of un- 
ravelling used for the calculations: as such, this phase is not a geometric object since it depends on 
non-physical parameters which are not related to the path followed by the density matrix during 
the evolution of the system. 
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I. INTRODUCTION 

The geometric phase is a property of a physical sys- 
tem which depends only on the path the system follows 
during its evolution, not on the details of the dynamics. 
Since the work of M. Berry Q], geometric phases have 
acquired a primary role in our understanding of many 
physical phenomena and have been subject to several 
experimental verifications The original idea, framed 
within the context of adiabatic and cyclic evolutions of 
isolated systems, has been generalized in various direc- 
tions |li@JajP, 81 ; of particular importance are all those 
efforts IE la S ISj iS] aiming at defining a geometric phase 
for open quantum systems, which are not mathematically 
described by a pure state {ipt) but in terms of a reduced 
density matrix pt undergoing a non-unitary evolution. 
Such proposals, besides been interesting on their own, 
could be important for the possible applications, e.g., to 
quantum computation IQJ. 

The mathematical framework within which the geo- 
metric phase for an open quantum system is defined is 
the following: 

1. One assumes that the effect of the environment 
on the quantum system is such that, under suitable 
approximations, the system can be effectively treated 
as an isolated system undergoing a non-unitary type of 
linear^ evolution: 



St : po St[/3o] = Pt, 



(1) 



which takes into account both the internal dynamics of 
the system and its interaction with the environment. Un- 
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^ As discussed in ref. ilU the evolution equation for the statistical 
operator must be linear, otherwise it can give origin to superlu- 
minal effects. 



der reasonable assumptions , which anyhow have been 
questioned in the literature (see e.g. JJJ), the map St 
can be taken of the quantum-dynamical-semigroup type, 
generated by the following class of equations: 



dp 
'dt 



N 

E 

n=l 



[L^LnP + pLlLn - 2LnpLl} ; 

(2) 

the self-adjoint operator H is usually identified with the 
standard Hamiltonian of the system, while the operators 
L„, together with the positive constant A, summarize the 
effect of the environment on the system. 

A general consequence of this type of approach is that 
a pure state \ipt) is usually mapped into a statistical 
mixture pt, so the problem arises of how to identify 
a geometric phase for the evolution of a density matrix pt . 

2. A common strategy which is used in the litera- 
ture for associating a geometric phase to the evolution of 
Pt is to formally map the density matrix into a statistical 
mixture of pure states \tp"), each of which is weighted 
with a probability p„ (t) : 



pt^m-),Pnm 



Pt 



= j2p„m^)m. (3) 



one can then use the standard definition of geometric 
phase for a pure state: 



It 



= Arg(V'o|^t)-Im / (^tMlV't), (4) 
Jo 

to associate a geometric phase also to pt ; this strategy as 
given fruitful results in the case of mixed states undergo- 
ing a unitary evolution ^ , while the case of non- unitary 
evolutions, in particular those associated to open quan- 
tum systems, is still under debate. In this second case, 
a tentative definition of geometric phase has been given 
via state purification 7] and the quantum-jump approach 
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A well-known characteristic property of relation (O 
is that the association between a density matrix pt and 
an ensemble {(|^"),p„(<))} is not one-to-one, but one- 
to-many T^]: in general, there are different ensembles, 
containing different vectors and different probabilities, 
which give rise to the same density matrix; moreover, 
such vectors and the corresponding probabilities evolve 
in completely different ways. In the state purification 
approach, this property of density matrices is related to 
the fact that there are different Kraus representations 
of the dynamical evolution of pt |l5j| . In the quantum- 
jump approach, the same property reflects the fact that 
there are different stochastic unravellings which end up 
to reproduce the same statistical operator [l^ . 

In this paper, we analyze the consequences of such a 
feature of density matrices for the definition of geometric 
phases, within the context of the stochastic unravelling 
formalism, which the quantum-jump approach belongs 
to. We will show that the approach followed in to 
identify a geometric phase strongly depends on the type 
of unravelling of pt'. as a consequence, such a phase is 
not a geometric object, because it depends also on non- 
physical parameters which are not related to the path 
followed by the density matrix during its evolution. 



II. STOCHASTIC UNRAVELLINGS 

A stochastic unravelling of a linear evolution St for a 
density matrix pt is defined as follows. Let us fix a prob- 
ability space (ri, P) and let us consider a stochastic 
evolution for statevectors: 



(5) 



which, for each different sample element w £ f2, asso- 
ciates a different statevector |^t(ci;)) to the same initial 
state l^/'o). One can then define the density matrix: 



Pt EE Y.p^MWt{^)){i^tH\], 



(6) 



where the symbol Ep denotes the average value with re- 
spect to the probability measure P; here above we have 
assumed that the initial state of the system is represented 
by a statistical mixture {(IV'o extra sum over 

the possible initial states appears in the definition of pt. 
The above relation defines a map: 



H ■ PQ 



(7) 



where pt is given by Eq. ^ . Now, if the map Et defined 
by Eqs. ((7J and JBJ coincides with the linear map Et, we 
say that Tt(w) is a stochastic unravelling of Sj. 

Among the other things, the above definition of 
stochastic unravelling implies that, when computing ob- 
servable quantities of a system which evolves according 
to the map Ej, e.g. the expectation value of a self-adjoint 
operator: {0)t = Tr[Opt], one can start with a stochas- 
tic unravelling (jSJ of St , then he computes the quantum 



expectation {i/jt{uj)\0\ipt{i^)) and finally he averages over 
the noise; this sequence of operations is legitimate since, 
by the definition of stochastic unravelling, one trivially 
has: 



Ep[(7At|0|Vt)] = Tr[OEt[|Vo)(^o|]]; 



(8) 



of course, if the initial state is a mixed state, an extra 
sum over the possible initial states, weighted with the 
corresponding probability distribution, has to be added 
at the left hand side of ©. 

In the literature, two types of stochastic unravellings 
have been proposed, one discrete and one continuous. 
The first one is the quantum jump approach which 
has been used in refs. Q to associate a geometric phase to 
the evolution of an open quantum system; the second one 
is given in terms of stochastic Schrodinger equations |l7j . 
The two approaches are similar and, in the following, we 
will resort to the second one, because it is more elegant 
from the mathematical point of view and easy to handle. 
We now briefly review it. 

The idea is simple: the stochastic evolution Tf (w) is as- 
sumed to be generated by a stochastic Schro ding er equa- 
tion, whose typical structure is the following [l7l llSl Il9j| : 



N 



Hdt + Xj2{Ln-rn,t)dWl' 



(9) 



where: 



A 

- V Y^^LiLn - 2L„r„,t + rlt) dt 



(10) 



with H and L„ defined as in (0); {n = 1,...,N) 
are N independent standard Wiener processes with re- 
spect to the measure P, which make Eq. a stochastic 
differential equation. 

Such kind of equations have been used in several con- 
texts: within the theory of quantum measurement, to de- 
scribe the effects of a repeated measurement on the evo- 
lution of a quantum system ; within collapse models, 
to provide a solution to the measurement problem [21^ : 
within the theory of open quantum system, as a math- 
ematical tool to efficiently simulate the evolution of an 
open system p^ . 

One of the fundamental properties 0| of Eq. is 
that the density matrix pt = Ep[\ipt) (iptl] solves Eq. lO, 
i.e. Eq. ^ represents a stochastic unravelling of the 
Lindblad-type equation jSJ. 

Note that Eq. lO is non-linear, but it preserves the 
norm of the statevector. There is a well-known way [T^ 
to linearize the equation, at the price of relinquishing the 
normalization condition; consider the following stochas- 
tic differential equation: 



N \2 ^ 

'-Hdt + XY,Ln - ^ E ^nLn dt 



(11) 
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the stochastic processes ^" are standard Wiener processes 
with respect to a new probabihty measure Q, whose re- 
lation to P will soon be established. 

The connection between the linear Eq. and the 
nonlinear Eq. is the following; given the solution 
of Eq. Hll|l for a initial condition {(po), if one performs 
the following two operations: 

1. Normalize the solution: \(j)t) — * \tpt) = \4>t) /\\\4't)\\, 

2. Make the substitution: 

Q Wr = ^r-2A frn,tds, (12) 

then the wavefunction Itpt) so defined is the solution of 
Eq. (|17|l for the same initial condition |V'o) = |0o)- More- 
over, one can further show that the two probability mea- 
sures P and Q are related as follows 

Ep[Xt] = EQ[{<j>t\<l>t)Xt], (13) 

where Xt is a stochastic process. 

A. Equivalent stochastic unravellings 
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FIG. 1: Interferometric scheme for measuring the total phase. 
|h} corresponds to the beam travelling in the horizontal di- 
rection, while |v) corresponds to the beam travelling in the 
vertical direction, x ^ variable phase shifter and B a mag- 
netic field. 

how the existence of different equivalent stochastic unrav- 
ellings affects the computation of the geometric phase. 
As an example, we now calculate the geometric phase as- 
sociated to the evolution of a spin particle in a constant 
magnetic field directed along the z-axis of a chosen ref- 
erence frame, while the spin is subject to dephasing. The 
quantum Hamiltonian is H = — /iScr^, and the effect of 
the environment is described by one Lindblad operator 
L = (7^; the corresponding Lindblad equation is {h = 1): 



Two stochastic unravellings T^^^ and t/^^ are said to 
be equivalent if they unravel the same evolution St. A 
very remarkable property is that there are infinite differ- 
ent but equivalent stochastic unravellings for practically 
all physically interesting Et; within the quantum jump 
approach, this issue is addressed, e.g., in |lj|; within the 
stochastic Schrodinger formalism, such a feature is less 
known, still very easy to show. As a matter of fact, sup- 
pose in Eq. Q we change the Lindblad operators i„ as 
follows: 



(14) 



where Cr, 



are arbitrary phase factors; clearly, 



Eq. ||2Jl does not change, while Eq. Q does change, since 
terms appear which are not proportional to L\Ln. Such 
a change is not as trivial as it may seem: as we shall 
see, stochastic equations with different values of (/?„ en- 
tail completely different evolutions for the statevector. 
Of course, there are other possible unravellings of Eq. 
besides those which can be obtained by a phase shift in 
the Lindblad operators, but for simplicity we consider 
here only these, since they are sufficient for the subse- 
quent analysis. 



d \^ 

—pt = ifiB [(7:,,pt] - — [(7^, Wz.Pt]] ■ (15) 



The initial spin state is taken equal to: 



IV'o) = COS-I+) + sin-|-), 



(16) 



where |+) and |— ) are the two eigenstates of Cz. Eq. 
becomes: 

d\ipt) = ifiBcr^dt + X{caz — cos If {az)t)dWt (17) 

1 

- y(cr^ - 2cCOS(y5 (CTz)tCTz +C0sV(0-z)t) IV't), 

with (crz)f = {ipt\'^z\''Pt)- In the above equation, we have 
included also the arbitrary phase factor c — e**^ which, 
as already discussed, does not appear in Eq. H15|l for 
the density matrix pt- We now compute the total and 
dynamical phases associated to the ensemble of vectors 
{\ipt) = |^/'t(w)},w e n} generated by Eq. ((T7|l . 



A. Total phase 



III. GEOMETRIC PHASE AND STOCHASTIC 
UNRAVELLINGS: AN EXAMPLE 

As already remarked, in refs. '9] the stochastic- 
unravelling approach has been used to associate a geo- 
metric phase to an open system; in this section we discuss 



To derive the correct formula for the total phase, we 
resort to an interferometric scheme like the one depicted 
in Fig. 1, i.e. a Mach-Zehnder interferometer with a vari- 
able phase shifter x in one of the two arms and the mag- 
netic field in the other. We call |h) and |v) the spatial 
part of the wavefunction when the beam travels along 
the horizontal and vertical direction, respectively; in the 
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subsequent analysis, we will neglect the contribution to 
the phase given by the free part of the evolution, since 
we assume that the length of the two arms of the inter- 
ferometer is the same. The initial state of the beam then 



is: 



®|h). (18) 



To find the evolution from the initial time t — to the 
final time t = when the beam comes out through 
the interferometer, we first consider the linear version 
of Eq. (^Tjl, which is^: 



d\<ft) = 



\2 



\^t). (19) 



By taking into account the effects of the two mirrors and 
partial beam splitters, one finds for the beam at time t^: 



+ e*^^^) e^^«*F- V(i+-^)*F cos ^ 1+) 



2 ' ' 



^ 2 



p-^(i+^ sin-l-) 
2 ' ' 



|v) 
(20) 



One can now compute the output intensity along |h): 

ItA^) = (*tp|[|h)(h|®Isp™]|«'tp) (21) 
= ^ + ^ I/«f(^)I cos(x + Arg(/tp(cj))), 

where \^t^) = |$tp)/|| |$tp) || and: 



g2ACtp cos if -itiBT ^^^2^ 



+ e 



sm — 
2 



(22) 



One can then identify the total phase, for each realization 
of the noise, as: 



7r('^)^Arg[/,(^)], 



(23) 



cannot be observed. The final outcome — i.e. the inter- 
ference pattern — consists of many repetitions of the ex- 
periment, accordingly the observable physical quantity is 
the average intensity = Ep[Jt(cL;)], which can be easily 
computed by taking into account relation (|13|l : 

It ^\ + \vT cos(x-7t'), (24) 

where: 



VT = \Ep[ftA^)]\ = 



is the output visibility, while: 
7$?' = ArgEp[/tp(w)] 



(25) 



Arg 



Arg EQ[(0o|0t)], 



(26) 



is the total phase difference; in particular, if B acts on 
the spin for a time T = ir/fiB, we have the standard 
result: 7^' = n. 

Note that the total phase 7.^' does not depend on (p, 
i.e. it does not depend on the specific choice of the un- 
ravelling used to make the calculations. Such a result is 
a consequence of property © of stochastic unravellings, 
i.e. of the fact that Eq. (|21|1 for the output intensity, 
when also the average over the noise is taken into ac- 
count, can be expressed as a function of the density ma- 
trix {It = Ep[It{uj)] = Tr[(|h)(h| (^IspiN)pt]), so that any 
dependence on the unravelling disappears. Such a result 
is then not a peculiar byproduct of the specific model 
taken into account, but a necessary mathematical con- 
sequence of the formalism. This fact can be seen in a 
different way: the total phase difference between the two 
arms of the interferometer is a physical quantity which 
can be experimentally measured; like all physical quanti- 
ties, it must be deducible from the master equation H15|l . 
so it does not have to depend on (p. 

As a final observation, we note that both the average 
visibility and the average total phase do not depend on 
A. This specific fact is a consequence of our simple model 
of open quantum system, according to which the noise is 
perfectly correlated among the two arms of the interfer- 
ometer. Of course such an assumption is not realistic, 
and it has been made only to simplify the calculations, 
since it does not affect the conclusion of our work. 



which depends not only on w, but also on i.e. on the 
type of unravelling of the master equation H15|l . Any- 
way, 7(°'(tj) as such is not a physical quantity because it 



^ In the following when we write |^t) and l^t), we mean that 
Eqs. 1171 and 1191 include not only the spin degree of freedom, 
but also the spatial one; when on the other hand we write |'(/>i) 
and \4>t), we mean that we are taking into account only the spin 
degree of freedom. 



B. Dynamical phase 

We now compute the dynamical phase 7^"" induced by 
the precession of the spin-system when interacting with 
the magnetic field; by using Ito calculus one finds 
from Eq. ATI) : 

(V'(i)MIV'W) = ifiB{a,)tdt + iXsmip{c7,)tdWt (27) 

\2 

— — [1 — {2c COS (fi — cos'^ip){(Tz) t] dt, 
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where only the imaginary part has to be taken into ac- 
count. For each reahzation of the noise, the dynamical 
phase is: 

ft rt 



{a z) s ds + X sin (f {az)sdWs 
Jo Jo 



+X'^ ain (p cos ip {az)lds, (28) 



which, like the total phase 7(°'(cj), depends not only on 
uj, but also on the unravelling of the master equation. 

We now compute the stochastic average 7^^" of 7^^"(tj), 
for which we need to know the statistical properties of 
both ((Tz)t and (crz)^: these can be easily computed by 
writing the corresponding stochastic differential equa- 
tions, both of which can be quite easily derived from 
Eq. l(T7|) . The equation for {az)t is: 

d{a,)t = 2Acos(^[l-(a.)2]dWt, (29) 

which tells us that since the Brownian increment dWt has 
zero mean, the average value of {<7z)t does not change in 
time: Ep[((Tz)t] = Ep[(cr2)o] = cos 6*. 

The stochastic differential equation for {az)t instead 

is: 



■Ep[(c7,)2] = 4A2cosVEp [1- (a,)?]' > 0, (31) 



d{azft = 4Xcosip{az)t[l~ {az)t]dWt 

+ AX'^cos'^ if [1- {a z)tf dt; (30) 

the first term on the r.h.s. does not contribute to the 
stochastic average, so one has: 

d 
dt 

which implies that 'Ep[{az)t] constantly increases in time, 
and in general (unless ip = 7r/2 + kn^k e Z) it stops 
increasing only when Ep [l — {o'z)t] — 0; when 
{az)t = Ij with the possible exception of a set of points 
w e il of measure 0. 

Concluding, the average dynamical phase 7^''^" = 
Ep[7^''"(w)], after a time T, is equal to: 

7^"" = /^Brcos6'-|-A^sin(/?cos(^ / Ep[{az)t]dt, (32) 

"'0 

which, contrary to what happens to the average total 
phase, still depends on ip, i.e. on the specific stochastic 
unravelling of the master equation. 

One could argue that, when computing the average 
dynamical phase, we should not average over the phase, 
but over the phase factor, i.e. we should compute 
Ep[exp(z7(''''°(ti;))] in place of Ep[7(''*'"(Li;)], and then ex- 
tract the argument; in this way we would take into ac- 
count the fact that a phase is defined modulus 2tt. The 
stochastic differential of exp{ijl°^ (ui)) is: 



i^B{(Jz)t + sin ip cos p{az)'t 



^ • 2 / \ 

-ysm p(az) 



Its average value cannot be explicitly computed, due to 
the dependence of both {az)t and {az)t on the noise; 
anyway, when taking the average, the dependence on ip 
in general does not disappear. E.g., if we take the triv- 
ial case in which the initial state is ji/'o) = \+), so that 
((jz)o = 1, then Eqs. and ||nU tell us that both {az)t 
and {crz)t remain equal to 1 for each realization of the 
noise; in such a case, the average value of exp(z7(°'(w)) 
at time T is: 



(34) 



Ep[exp(i7^°*(cj))] =e-- 
and its argument clearly depends on (p. 

C. Geometric phase 



The geometric phase 7(°° is the difference between the 
total and the dynamical phase. For each realization of the 
stochastic process Wt, one has from Eqs. and (^5)) '^: 



7-°(u;) = Arg[/t(u;)] -/iS / {cTz),ds (35) 

Jo 

- Xsinp {cFz)s dWs 



— X svnpcosip I {az)gds 
Jo 

(/((cj) is defined in H22|l ). which clearly depends on the 
type of unravelling. Its average value is: 



7f° = ArgEp[/t(c^)] 



Ep[7^^"(cj)] or 
Arg[Ep[exp(z7^^"(o 



(36) 



Also the average geometric phase depends on the type of 
stochastic unravelling of the master equation, whichever 
way the average dynamical phase is computed. This is 
the main result of our paper. 



IV. DISCUSSION AND CONCLUSIONS 

The geometric phase of an open quantum system 
should be a quantity depending only on the path fol- 
lowed by the density matrix pt in its state space; we 
have seen that such a phase, when computed by means 
of stochastic unravellings — as done in ref. Q — depends 
on the type of unravelling, both for single realizations of 
the noise (Eq. H35|ll and for its average value (Eq. H36(l ). 
This fact as two important consequences: 



+ [iXsinp{(7z)t] e*'^'°'(")dWt. 



(33) 



^ In Eqs. I35i and 1361 . t refers to the time during which the beam 
travels through the interferometer, while T is the time during 
which the spin interacts with the magnetic field. 
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1. First of all, the phase defined in 1)35(1 and (|36|l is 
not a geometric object, since it depends also on ip 
which has nothing to do with the path followed by 
Pt during its evolution. 

2. Worse than this, such a phase is not even an object 
somehow related to a physical quantity, because ip 
itself has no physical meaning since it only selects 
one of the infinitely many equivalent stochastic un- 
ravellings which can be used. 

The conclusion is that the stochastic-unravelling method 
does not lead to a sensible definition of geometric phase. 

One could say that different stochastic unravellings 
might correspond to different ways to perform the mea- 
surement or to monitor the environment |16| and then 
that different values for the phase correspond to differ- 
ent ways to measure the system; anyway, this is not here 
the case: in our example we have taken a standard in- 
terferometer where the output intensity is measured in a 
standard and unique way. Nevertheless, different unrav- 
ellings can still be taken into account. 

Note that the dependence on ip comes only from the 
dynamical phase, not from the total phase. As already re- 
marked, this is not a consequence of the specific model of 
open quantum system we are considering here, but a di- 
rect consequence of the fact that a total-phase difference 
is a measurable quantity and as such must be deducible 
from the density matrix pt, which docs not depend on ip. 
On the contrary, the dynamical phase and thus also the 
geometrical phase, is not directly observable so — at least 
from the mathematical point of view — it can depend on 

as it happens here. 

What is the mathematical origin of the dependence of 
^l°°(uj) and on (pi Its unravelling-dependence does 
not come from the total phase but from the dynamical 
component, and for the following reason: by definition, 
^'^^"{uj) is not a function of li/it) at the considered time, 
but a function of the whole history of IV's), from s — io 
s = t, i.e. it depends on the whole trajectory followed by 
the statevector. Now, it is easy to see that for different 
unravellings the trajectories followed by the statevector 
are radically different. For example, when ip = Tr/2 + 
kiT with fc € Z, Eq. (|29|) tells that {a^) is constant in 
time, for each realization of the stochastic process: this 
implies that the projection of the spin vector along the 
magnetic field does not change in time, i.e. the vector 
rotates always along the same circle on the Bloch sphere. 
On the other hand, when (/? 7^ 7r/2 + /ctt with fc G Z then, 
as we have already discussed in connection with Eq. (|31|l . 
(ctz)^ approaches the value 1 for t — > 00, i.e. the variance 
Var[cr2] = {af) - (cr^)^ = 1 - (cr^)^ of the operator CTz 
approaches zero: this means that the statevector is driven 
towards one of the two eigenstates of CTz, thus changing 



the projection of the spin vector along the magnetic field 
|23l |. As a consequence, being the trajectories followed 
by the statevector so strongly dependent on the kind of 
unravelling, there is no need for the dynamical phase to 
be unravelling independent, as it actually occurs. 

A different way to see what happens is the following: 
the relation {ipt\d\ipt) — defines the parallel transport 
condition, and for different unravellings one has different 
inequivalent parallel transport conditions, thus different 
definitions of a geometric phase. 

Another interesting question is about the physical rea- 
son for such a dependence of the geometric phase on 
ip. Someone has argued that the problem arises be- 
cause the master equation used to model the effect of 
the environment is of the Lindblad type. Since the Lind- 
blad equation is only an effective equation approximating 
an otherwise too complex system, its validity is limited 
and it could not be suitable for computing the geometric 
phase. We think that this is not the case: the source of 
all troubles derives from the fact that there are different 
equivalent stochastic unravellings associated to the same 
evolution Et which determines different evolutions for the 
statevector, thus different parallel transport conditions; 
such a feature is not an exclusive property of the Lind- 
blad equation (it is not even a mathematical consequence 
of it), but has a more general character. 

To summarize, the stochastic-unravelling approach 
used in does not produce a phase which is geomet- 
ric, i.e. which depends only on the trajectory followed 
by Pt during the evolution; it depends also on the spe- 
cific choice of the unravelling used for the calculations, 
which by itself has no particular physical meaning. This 
difficulty can be in principle overcome by fixing the un- 
ravelling to be used for computing the geometric phase, 
as often implicitly done in the literature, but this pro- 
cedure cannot be satisfactory for two reasons: first, it 
obviously does not remove the fact that the definition is 
mathematically unravelling-dependent; second, there is 
no fundamental physical reason why to choose one un- 
ravelling in place of another, since they are all on the 
same footing. 

In ref. it has been stated that, within the state pu- 
rification approach of '7| , different Kraus representations 
may lead to different values for the geometric phase: if 
so, then our criticism apply also to the approach of 



V. ACKNOWLEDGEMENTS 

The work of A.B. was supported by the Marie Curie 
Fellowship MEIF-CT-2003-500543. The work of E. I. was 
supported by the "Consorzio per la Fisica — Trieste" . 



[1] M.V. Berry, Proc. R. Soc. Lond. A 392, 45 (1984). 



[2] Geometric Phases in Ptiysics, edited by A. Shapere and 
F. Wilczek (World Scientific, Singapore, 1989). 



7 



[3] R. Tycko, Phys. Rev. Lett. 58, 2281 (1987). R. Bhandari 
and J. Samuel, Phys. Rev. Lett. 60, 1211 (1988). R.Y. 
Chiao, A. Antaramian, M.K. Ganga, H. Jiao and S.R. 
Wilkinson, Phys. Rev. Lett. 60, 1214 (1988). D. Suter, 
K.T. Mueller and A. Pines, Phys. Rev. Lett. 60, 1218 
(1988). 

[4] Y. Aharonov and J. Anandan, Phys. Rev. Lett. 58, 1593 
(1987). J. Samuel and R. Bhandari, Phys. Rev. Lett. 60, 
2339 (1988). A.K. Pari, Phys. Rev. A 52, 2576 (1995). 

[5] A. Uhlmann, Rep. Math. Phys. 24, 229 (1986). D. Gam- 
liel and J.H. Freed, Phys. Rev. A 39, 3238 (1989). D. 
ElUnas, S.M. Barnett and M.A. Dupertuis, Phys. Rev. A 
39, 3228 (1989). K.M. Fonseca Romero, A.C. Aguir and 
M.T. Thomaz, Physica A 307, 142 (2002). 

[6] E. Sjoqvist, A.K. Rati, A. Ekcrt, ,J.S. Anandan, M. Eric- 
sson, D.K.L. Oi and V. Vcdral, Phys. Rev. Lett. 85, 2845 
(2000). 

[7] M. Ericsson, E. Sjoqvist, J. Brannhmd, D.K.L. Oi and 
A.K. Rati, Phys. Rev. A 67, 020101(R) (2003). D. Kult 
and E. Sjoqvist, Int. J. Quantum Comp. 2, 247 (2004). 
K.-P. Marzlin, S. Ghose, and B.C. Sanders, Phys. Rev. 
Lett. 93, 260402 (2004). 
[8] D.M. Tong, E. Sjoqvist, L.C. Kwek, and C.H. Oh, Phys. 

Rev. Lett. 93, 080405 (2004). 
[9] A. Carollo, L Fuentes-Guridi, M.F. Santos and V. Ve- 
dral, Phys. Rev. Lett. 90, 160402 (2003); Phys. Rev. 
Lett. 92, 020402 (2004). X. X. Yi, D. P. Liu, and W. 
Wang, New J. Phys. 7, 222 (2005). 

[10] P. Zanardi and M. Rasetti, Phys. Lett. A 264, 94 (1999). 
J. A. Jones, V. Vedral, A. Ekert and G. Castagnoli, Na- 
ture 403, 869 (2000). J. Pachos, P. Zanardi and M. 
Rasetti, Phys. Rev. A 61, 010305(R) (2000). A. Ekert, 
M. Ericsson, P. Hayden, H. Inamori, J. A. Jones, D.K.L. 
Oi and V. Vedral, J. Mod. Opt. 47, 2051 (2000). G. Falci, 
R. Fazio, G.M. Palma, J. Siewert and V. Vedral, Nature 
407, 355 (2000). L. M. Duan, J.L Cirac, and P. ZoUer, 
Science 292, 1695 (2001). 

[11] N. Gisin, Phys. Lett. A 143, 1 (1990). 

[12] G. Lindblad, Rep. Math. Phys. 10, 393 (1976). 

[13] D. Kohen, C. C. Marston, and D. J. Tarmor, 
J. Chem. Phys. 107(13), 5236 (1997). 

[14] M.A. Nielsen and LL. Chuang, Quantum Computation 
and Quantum Information, Cambridge University Press, 
Cambridge (2000). 

[15] K. Kraus, States, Effects and Operations, Springer- 



Vcriag, Berlin (1983). 

[16] H. Carmichael, An Open Systems Approach to Quantum 
Optics, Springer- Verlag, Berlin (1993). 

[17] A. Barchielli, in Contributions in Probability - In mem- 
ory of Alberto Frigerio, edited by C. Cecchini, (Forum, 
Udine, 1996). 

[18] L. Arnold, Stochastic Differential Equations: Theory and 

Applications, Wiley, Now York (1974). 
[19] G.W. Gardiner, Handbook of Stochastic Methods, 2nd 

Ed., Springer, Berlin (1997). 
[20] V.P. Belavkin, Journ. Math. Phys 31, 2930 (1990). A. 

Barchielli and V.P. Belavkin, Jounr. Phys. A 24, 1495 

(1991). V.P. Belavkin and P. Staszewski, Phys. Rev. A 

45, 1347 (1992). 
[21] P. Pearle, Phys. Rev. D 13, 857 (1976); Phys. Rev. D 

29, 235 (1984); Phys. Rev. Lett. 53, 1775 (1984). L. 

Diosi, Phys. Lett. A 132, 233 (1988); Journ. Phys. A 

21, 2885 (1988); Phys. Rev. A 40, 1165 (1989); Phys. 

Rev. A 42, 5086 (1990). G.C. Ghirardi, P. Pearle and 

A. Rimini, Phys. Rev. A 42, 78 (1990). A. Bassi and 

G. C. Ghirardi, Phys. Ropt. 379, 257 (2003). L. Diosi, 
Phys. Lett. A 132, 233 (1988); Journ. Phys. A 21, 2885 
(1988); Phys. Rev. A 40, 1165 (1989); Phys. Rev. A 42, 
5086 (1990). S.L. Adler and L.P. Howritz, Journ. Math. 
Phys. 41, 2485 (2000). S.L. Adler and T.A. Brun, Journ. 
Phys. A 34, 4797 (2001). S.L. Adler, Journ. Phys. A 35, 
841 (2002) S.L. Adler, Quantum Theory as an emergent 
phenomenon, Cambridge Univorsityn Press, Cambridge 
(2004). 

[22] N. Gisin, Phys. Rev. Lett. 52, 1657 (1984); Helvet. Phys. 
Acta 62, 363 (1989). N. Gisin and L Percival, Journ. 
Phys. A 25, 5677 (1992); Journ. Phys. A 26, 2233 (1993). 
N. Gisin and M. Rigo, Journ. Phys. A 28, 7375 (1995). 

H. P. Brouer and F. Pctruccionc, Journ. Phys. A 31, 33 

(1998) . H.P. Breuer and F. Petruccione, in Open Systems 
and Measurement in Relativistic Quantum Theory, H.P. 
Breuer and F. Petruccione eds., Springer- Verlag, Berlin 

(1999) . H.P. Breuer and F. Petruccione, in Chance in 
Physics: Foundations and Perspectives, J. Bricmont et 
al. eds. Springer- Verlag, Berlin (2001). H.P. Breuer and 
F. Petruccione, The Theory of Open Quantum Systems, 
Oxford University Press, Oxford (2003). 

[23] A. Bassi, and E. Ippoliti, Phys.Rev. A 69, 012105 (2004). 
[24] Comment by one of the referees. 



